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■ The precise observation of the angle-frequency spectrum of light filaments in water reveals a 
/^S ' scenario incompatible with current models of conical emission (CE). Its description in terms of 

linear X-wave modes leads us to understand filamentation dynamics requiring a phase- and group- 
matched, Kerr-driven four- wave-mixing process that involves two highly localized pumps and two 
X-waves. CE and temporal splitting arise naturally as two manifestations of this process. 
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PACS numbers: 190.5940, 320.2250 

Filamentation of intense light pulses in nonlinear me- 
dia has attracted much interest ever since first exper- 
imental evidences in the early Ws (Q and references 
therein) . Owing to the very high intensities reached dur- 
ing the process, several nonlinear phenomena, e.g., multi- 
photon absorption, plasma formation, saturable nonlin- 
ear response, stimulated Raman scattering etc., occur 
in addition to the optical Kerr effect. Indeed, the fila- 
ment regime is enriched by peculiar phenomena like pulse 
splitting, self-steepening, shock-wave formation, super- 
continuum generation, and conical emission (CE) 0. In 
media with normal group velocity dispersion (GVD), no 
matter if of solid, liquid or gaseous nature, CE accom- 
panies filamentation, producing radiation at angles that 
increase with increasing detuning from the carrier fre- 
quency 0, Q. In spite of the generality of the process, 
a clear understanding of the interplay between CE and 
filament dynamics is still missing. Only recently, Kolesik 
et al. have proposed an interpretation of filamentation 
dynamics in water on the basis of pulse splitting and dy- 
namic nonlinear X waves at the far field |5j, in which 
the double X-like structure observed in simulated angle- 
frequency spectra arises from the scattering of an incident 
field at the two main peaks of the split material response 
wave. 

Originally CE in light filaments was interpreted in 
terms of the modulation instability (MI) angle-frequency 
gain pattern of the plane and monochromatic (PM) 
modes of the nonlinear Schrodinger equation (NSE) |y,l3- 
Measurements at large angles and detunings from the 
carrier frequency gave in fact results fairly compatible 
with this interpretation H, ■ In the present work, ow- 
ing to the use of a novel imaging spectrograph technique 
|3j, we have been able to observe for the first time the 
CE in the region of small angles and detunings. The re- 
sults clearly indicate a scenario not compatible with the 



MI analysis of PM modes. Our description by means of 
the spectra of the stationary linear X-waves supported 
by the medium, indicates that the strong localization of 
the self-focused field plays a crucial role in the substantial 
modification experienced by the MI pattern. We propose 
a simple picture in which the latter results from the para- 
metric amplification of two weak X-waves by the strong, 
highly localized pump. Supporting this interpretation, 
we are able to derive, from the matching condition among 
the interacting waves, a simple analytical expression [Eq. 
(0}] that accurately determines the overall CE structure, 
and predicts also the splitting velocity of the pump wave 
as a function of the peak intensity reached at collapse. 

The experiment was carried out using a 3 cm long cell 
filled with water as a bulk nonlinear Kerr medium. Beam 
filamentation was induced by launching 200 fs pulses 
at carrier wave length 527 nm, delivered from a mode- 
locked, chirped-pulse regeneratively amplified, Nd:glass 
laser with a 10 Hz repetition rate (Twinkle; Light Con- 
version Ltd.). The pulse energy was controlled by use of 
a half-wave plate and a Glan- Taylor polarizer. A spatial 
filter guarantees high beam quality before focusing with 
a / = 50 cm lens placed at 48 cm from the cell input 
facet, with a beam diameter at half-maximum equal to 
100 /im. The output facet is then imaged with a 4-/ 
telescope onto the rear focal plane of the lens (fp = 15 
cm) used to obtain the spatial Fourier transform of the 
filament generated in the Kerr sample. The input slit 
of the imaging spectrograph (MS260i, Lot-Oriel) with 
a high resolution 1200 lines/mm diffraction grating is 
placed at a distance fp from the Fourier lens so as to 
reconstruct the angle-wavelength (0-X) spectrum of the 
filament on a commercial CCD 8- bit camera (Canon), 
placed at the output plane of the spectrograph. With 
the change A = 2ttc/(ojq + ft), the angle- frequency (6- 
Q) spectrum can also be obtained, Q being the detuning 
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FIG. 1: Measured 0-X spectrum for E ln = 2 uj. 

from the carrier frequency uiq. More details of the exper- 
imental layout may be found elsewhere Q ■ We recorded 
only single-shot spectra in order to avoid averaging ef- 
fects due to possible shot-to-shot fluctuations in the input 
pulse energy. Moreover we highly saturated the central 
peak of the spectrum so as to highlight the less visible 
surrounding structure. For input energies E in < 1.8 fiJ, 
no CE, or clear X-like features were observed. Figure 
shows an example of 9-X spectrum at E lVL = 2 /iJ. The 
CE appears as distinctly separated red- and blue-shifted 
X-shaped tails. This pattern remains very similar with 
increasing input energy up to E- ln ~ 4 /it J, while further 
increase produces a slowly deteriorating picture, with a 
modulated intensity pattern that extends to nearly all 
recorded values of 9 and A. We interpret this deteriora- 
tion as due to the onset of local breakdown in the water 
sample. 

CE emission is commonly accepted to arise from the 
interplay of diffraction, dispersion and nonlinear material 
response, the simplest model that accounts for it being 
the cubic NSE with normal GVD [10 

f^J-V^-f + (1) 
oz 2ko 2 ar z c 

Here, A(x,y,T, z) is the complex envelope of the wave 
packet E — Aexp(ikoz — iuj^t) of carrier frequency 
ujq, = d 2 /dx 2 + d 2 /dy 2 is the two-dimensional 

Laplace operator perpendicular to the propagation di- 
rection z, t = t — k' Q z is the local time, c the speed of 
light in vacuum, 712 the nonlinear refractive index, and 
fcg = [d m k(uj)/duj m ]uj , where is k(cu) the frequency- 
dependent propagation constant. Mathematically, the 
9-fl spectrum where CE is observable, is directly related 
to the Fourier spectrum of the envelope. If for instance, 
A(r, t, z) [r = (x 2 + J/ 2 ) 1//2 ] is a cylindrical symmetric 
complex envelope, and A(K±, O, z) [K± = (k^+k 2 ) 1 ^ 2 ] is 
its spatiotemporal Fourier transform, then the 9-fl spec- 
trum is given by A(k$9, f2, z), where 9 = K±/ko is the 
propagation angle with respect to the z axis. 



An accepted approach for understanding the structure 
of the CE relies upon the evaluation of the MI gain pro- 
file of the PM modes of the NSE 0] . In Kerr self- focusing 
media with normal dispersion, the perturbations to a PM 
mode that grow at maximum rate are those whose spa- 
tiotemporal frequencies (K±, Q) are related by 

K X (Q) = ^/fc fco fi2 + 2fc oA > 0), (2) 

i.e., lie on a hyperbola on the K±-£l plane featuring an 
angle gap (K± gap) 0. In this respect, Luther et al. 
proposed an intuitive picture that assumes the largest 
MI gain to occur at angles and frequencies fulfilling the 
linear phase-matching condition of the four-wave mixing 
(FWM) process supported by the Kerr response • Un- 
der this hypothesis, the asymptotic linear approximation 
K±($i) ~ ko fcp I f2 1 was re-obtained, and the observed 
discrepancies were attributed to the nonlinear phase shift 
produced by the PM mode on the weak perturbation. 

Precedingexperimental observations of CE in borosil- 
icate glass Q and ethylene glycol were indeed inter- 
preted to present a hyperbolic structure with an angle 
gap, which was also attributed to the nonlinear phase 
shift. This trend is also visible in our measurements 
(Fig. ^) as an intersection of the two X arms at a 
non-zero angle. If we attempt to fit Eq. J2J (taking 
f3 and fcg as free parameters) to the experimental data 
[Fig'EI a ), dotted line and closed circles, respectively] we 
obtain fcg = 0.053 ± 0.03 fs 2 /^m, which firmly sup- 
ports the predictions regarding the asymptotic slope 
(fco = 0.055 fs 2 //xm at 527 nm however, the angle- 

gap fitting strongly departs from experimental data in 
the non-asymptotic region, which strongly calls for a fit- 
ting to a hyperbola with opposite curvature, that is, to 
a hyperbola with a frequency gap. This observation ad- 
vocates a novel interpretation of the origin of CE. 

To understand our proposal, note first that the usual 
MI interpretation of the seemingly angle-gap hyperbolic 
CE can be equivalently reread in terms of the excitation 
of a weak, linear X-wave mode by a strong, nonlinear 
PM pump. Indeed, Eq. J5J) represents also the K±-fl 
spectrum of an X-wave mode of the medium [Tlj . that is, 
a diffraction- and dispersion-free pulsed Bessel beam in 
which cone-angle-induced dispersion [6(d) ~ K±(n)/k ] 
and material GVD balance mutually. The spectrum in 
Eq. (J2J) belongs to a wave-mode of carrier frequency Qq 
and shortened axial wave vector feo — /3. In this reading, 
the excitation of a X-wave mode with angle-gap (j3 > 0) 
is a consequence of the nonlinear phase matching be- 
tween PM pump and X-wave for maximum efficiency of 
the interaction: taking into account the nonlinear cor- 
rections to the wave vectors of both pump and X-wave, 
phase matching imposes fco + Afc = ko — (3 + 2Ak, where 
Afc = ujQn2l/c is the positive nonlinear phase shift (in 
self-focusing Kerr medium) for the PM pump of intensity 
/, and where 2Ak is the corresponding nonlinear correc- 
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tion to the weak perturbation 0, , leading immediately 
to ft = UMniljc > (as predicted by the standard MI 
analysis p). 

We also note that a medium with GVD [k(Cl) = k + 
fc f2 + fc '£l 2 /2] supports more general X-wave modes with 
strengthened and shortened wave vectors (/3 negative and 
positive), as well as shifted carrier frequency u>n + th eir 
most general K±-H, spectrum being expressed by 



(/3 = 0). 



(3) 



Their phase and group velocities are given by = 
(lu + Cl)/[k(Q) - (3] and v^) = l/k'(Q), respectively, 
which can take both subluminal and superluminal val- 
ues. Equation @ describes a two-parameter family of 
hyperbolas, sharing the same asymptotic slope y'fcofc^, 
but admitting angle gaps {(3 > 0) and frequency gaps 
(3 < 0), as well as positive and negative frequency shifts 
f2. X-wave spectra appear then as a suitable tool for the 
description of the CE features. 

The apparent frequency-gap hyperbolic form of the ob- 
served CE, and the actual absence of such a gap, suggest 
that the CE cannot be described in terms of a single X- 
wave. This leads us to interpret the two observed X tails 
as belonging to two different X waves both featuring fre- 
quency gaps, and proceed to fit Eq. J3J with two indepen- 
dent sets of parameters f2 and (3 (and k' ' = 0.056 fs 2 /^m) 
to the experimental data. This procedure also finds moti- 
vation in the X-X structure of the numerically evaluated 
angular spectra of light filaments in water, which was 
interpreted to be a consequence of pulse splitting on X- 
wave generation Fittings with Vl = +0.33, —0.33 fs _1 
and (3 — —2.2, -2.2 mm" 1 (dashed and continuous lines, 
respectively) reproduce accurately the structure of the 
CE [Fig. Ufa)]. The interpretation in term of X- waves 
leads moreover to suspect the existence of two additional, 
not yet observed X arms in the K-VL spectrum at wave- 
lengths 450 and 635 nm [Fig. E^a)]. 

To verify this prediction, we performed additional mea- 
surements with a lower resolution, 300 lines/mm diffrac- 
tion grating in order to cover a wider spectral region, as 
shown in Fig. Efb). Alongside the central tails (around 
527 nm), we clearly observe the generation of new fre- 
quencies in the vicinity of 450 and and 635 nm, as ex- 
pected. The reduced visibility of the red tail at 635 nm 
can be attributed to the stronger linear absorption at this 
wavelength. In fact, we had to increase the input pulse 
energy up to E in = 3 jit J in order to enhance the overall 
intensity of the new tails. 

The question then arises of the mechanism responsi- 
ble for the described spectral structure. Here we pro- 
pose that the double X-like CE is a result of a nonlinear 
FWM interaction — the lowest order process supported 
by the nonlinear Kerr response — , in which two signal 
and idler X-waves experience a parametric amplification 
by a strong, highly localized (i.e., non PM) pump. If as 
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FIG. 2: (Color online) a) Closed circles: 6-Q. distribution 
of peal fluence in Fig. dotted line: best fit of the entire 
set of experimental data with Eq. (J2J; red-dashed line: best 
fit of the Q. < data with with Eq. 10; blue-full line: best 
fit of the Q, > data with with Eq. (J^J b) Experimentally 
measured 0-X spectrum for Ei n = 3 (ij, showing the Stokes 
and anti-Stokes sidebands. 



consistently with the self-focusing dynamics, we consider 
a strongly localized pump rather than a PM wave, it is 
expected that nonlinear phase modulation does not in- 
fluence significantly an extended object such an X-wave 
mode, whose axial wave vector is solely determined by 
its cone angle, that is, by its weak, longstanding tails far 
from the pump. In fact, in the X-wave the energy does 
not flow along the axial direction, but along a conical 
surface, which prevents pump-induced cross-phase accu- 
mulation to take place. 

Under this assumption, phase matching between a lo- 
calized pump and a single X-wave at same frequency 
ujq would require fco + Afc = ko — (3, a condition that 
contrarily to the case of a PM pump, is satisfied by a 
frequency-gap X-wave mode (f3 — —Ak < 0). More gen- 
erally, if k = ko + Afc is the nonlinear pump wave num- 
ber, and k s ^i — k(fl St i) — (3 S ^ are two, signal and idler, X 
wave numbers, where k(fl s ^i) = fco + k fl St i + fc 'S7 2 i /2, 
and are their carrier frequency shifts, the conditions 
of energy and momentum conservation, fi s = — and 
2fc = k s + ki, for maximum efficiency of a FWM pro- 
cess involving two (degenerate) highly localized pumps 
and the signal and idler X waves, leads to the relation 
$s + A — fed'fifi = — 2Afc. Among all possible couples of 
X-wave modes satisfying this relation, those whose spec- 
tra cross the pump, located around £1 = 0, K± = 0, are 
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FIG. 3: Left: Values of f3 s ,i obtained from X-wave fitting to 
the numerical K±— fi spectra for increasing pump intensities 
I (open squares), and prediction of Eq. gj (solid curve). 
Right: GVM between the splitting pulses obtained from direct 
space-time numerical simulation (open circles), and between 
the X- waves as given by 2\J k^u^-nil j c (solid curve). 



the most energetically favored, since these X-waves will 
not need to grow from noise, but from the more effec- 
tive pump self-phase modulation. This condition was 
also found in 12] to result in the largest gain in the 
case of a single X wave excited by a travelling pump. 
The X-wave-pump crossing condition leads, from Eq. 
©, to (3 St i = — fco'Ogi/2, and then, on account that 
Afc = uJori2l/c, to 



-Wo^2^/2c . 



(4) 



Accordingly, the two X-waves will present symmetrically 
shifted carrier frequencies f^. s = io§niI j ck§, and 
frequency-gaps of width 2|f2j iS |, from fi = towards the 
Stokes and anti-Stokes bands, as observed in the experi- 
ment. 

Eq. (@J predicts also a precise dependence of the whole 
CE structure on the pump intensity J, whose validity we 
have tested numerically. We solved the NSE with 
fc ' = 0.056 fs 2 /^m, n = 1.33 and n 2 = 1.6 x 10~ 16 
cm 2 /W, for an input 200 fs-long, 100 ^m-wide (FWHM) 
Gaussian wave packet. Since CE is seen to develop explo- 
sively at the collapse, we identified the pump intensity / 
with the absolute peak intensity reached during propaga- 
tion. In order to attain different values of this intensity 
without changing the rest of parameters, we added to the 
second member of the NSE (JJJ the nonlinear dissipation 
term —P^ K '\A\ K A/2, with K = 3 (three-photon absorp- 
tion at 527 nm in water), and with (}( K ) ranging from 
1.2 x 10~ 23 to 8 x 10~ 25 cm 3 /W 2 . The K x - spec- 
tra obtained numerically are symmetric with respect to 
£1 = due to the approximations involved in the NSE, 
and present a nearly invariant X-X structure beyond the 
collapse point (in spite of the quickly decreasing peak 
intensity /), with one tail of each X wave nearly pass- 
ing through £1 = 0, as observed experimentally. The 



spectrum at 3 cm from the input plane was then easily 
fitted with two, signal and idler, material X-waves modes 
crossing il = [i.e., with fc ' = 0.056 fs 2 //mi, and with 
f2 Sl i = ±(— 2/3 Si j/fco) 1 ' 2 ], fis,i being then the only one free 
parameter. Figure shows the values of /3 Si i, obtained 
from the best fits to the spectra, versus pump intensity 
/ (open squares). The excellent agreement with the pre- 
dicted P s ,i~I dependence in Eq. Q (solid line) strongly 
supports the FWM analysis. 

Extending further our interpretation, we conjecture 
that the well-known phenomenon of pulse temporal split- 
ting in filamentation with normal GVD, usually de- 
scribed from a purely temporal perspective, is a conse- 
quence of higher-order matching among the interacting 
waves — in our case, of group matching. Indeed, the zero- 
order, phase-matching condition Q entails that the two 
X-waves must travel at different group velocities v^} = 
l/k'{yi S: i) — l/[fcg -I- fcgf2 Si i], and therefore split apart in 
time one from another with a group velocity mismatch 
(GVM) l/vi 9) - l/v[ 9) = 2fc '|Q s>i | = 2y/k%u> n2l/c pro- 
portional to the square root of the pump intensity. Group 
matching among the interacting waves is then better at- 
tained if the two pump waves, breaking their initial de- 
generacy, split also to co-propagate with the X-waves. In 
this view, pulse splitting is not a mere collapse-arresting 
mechanism Q, which further determines, taken it for 
granted, the generation of two X-waves || ; instead, pulse 
splitting emerges here as a particular feature of the phase 
and group matched wave configuration most favored by 
the FWM nonlinear interaction inherent to the Kerr NSE 
dynamics. In other words, X-wave instability or paramet- 
ric amplification, and splitting instability are not inde- 
pendent phenomena, but two aspects of the same process. 
To sustain this hypothesis, we have evaluated the GVM 
of the split pumps in the same numerical simulations as in 
the preceding paragraph. To do so, we plotted the tem- 
poral profiles at different propagation distances beyond 
collapse, and obtained the GVM as the slope of a linear 
fitting to the delay between the two temporal peaks ver- 
sus propagation distance. Figure [3] shows that the pump 
GVM (open circles) does not depart more than 10 per- 
cent from the predicted X-wave GVM (solid curve), and 
that both follow a similar dependence with the square 
root of the intensity / at collapse. 

In summary, our experiments demonstrate that CE 
emission in the angular spectrum of filaments is not in- 
terpretable in the frame of MI of PM modes of the NSE. 
Strong localization of the self- focusing pulse substantially 
modifies the MI pattern, which finds accurate description 
in terms of linear X-wave modes of the medium, and sim- 
ple explanation as a result of a dominant, phase-matched 
FWM mixing process supported by the NSE dynamics 
between two highly localized, strong pump waves and 
two amplifying weak X-waves. Pulse temporal splitting 
emerges in this model as the necessary temporal dynam- 
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ics for preserving group matching among the interacting 
waves. 
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